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Abstract 

The soft-gluon resummation exponents G N in moment space are investigated for the quark 
coefficient functions in deep-inelastic structure functions and the quark- ant iquark contri- 
bution to the Drell-Yan cross section da/dM. Employing results from two- and three-loop 
calculations we obtain the next-to- next-to-leading logarithmic terms a s (a s In N) n of G N 
to all orders in the strong coupling constant a s . These new contributions facilitate a 
reliable assessment of the numerical effect and the stability of the large- N expansion. 



Deep-inelastic lepton-hadron scattering (DIS) and Drell-Yan (DY) lepton pair production 
in hadronic collisions are among the processes best suited for probing the short-distance 
structure of hadrons. The subprocess cross sections (coefficient functions) for these pro- 
cesses in perturbative QCD receive large logarithmic corrections, originating from soft- 
gluon radiation, at large values of the scaling variables x corresponding to large Mellin 
moments N. These corrections have been resummed [l], [| to all orders in the strong 
coupling constant a s up to next-to-leading logarithmic accuracy (technically defined in 
eq. fllOD below). As the next-to- leading contributions are often hardly suppressed against 
the leading terms, it is important to extend the resummation to the next-to-next-to- 
leading logarithms. In this paper we present corresponding results for the DIS structure 
functions Fi^^x, Q 2 ) (where Q 2 represents the resolution scale) and for the Drell-Yan 
cross section da/dQ 2 (where Q stands for the invariant mass M of the lepton pair). 

The soft-gluon resummation of the DIS (P = 1) and DY (P = 2) quark coefficient 
functions Cp(Q 2 ) in iV-space is given by |2| 

Cp(Q 2 ) = 9p,o(Q 2 ) ■ ex P [Gp{Q 2 )} + 0{N-Hn n N) . (1) 

The functions C N , g Q and G N also depend on the factorization scale fif and the renormal- 
ization scale fi r , a dependence which we will often suppress for brevity. For the quantities 
under consideration the contributions go collecting the iV-independent terms are known 
up to second order in as from the two-loop calculations performed in refs. ||. In the 
standard MS renormalization and factorization scheme employed throughout this paper 
the exponents Gp in eq. (|1|) can be written as 

G?(Q 2 ) = lnA ? (Q 2 ,/4) + In J q (Q 2 ) + In Af(Q 2 ) , 

G%{Q 2 ) = 2 1nA 9 (Q 2 ,/4) + lnA:f(Q 2 ) . (2) 

Closely following the notations of ref. ||, the components entering eq. (|2]) are 

rl z N-l _ i r (l~z) 2 Q 2 AfP- 

lnA,(Q 2 ,^) = / dz— / -\A(a s (q 2 )) (3) 

Jo i — z Jfii q 

collecting of effects of soft-gluon radiation collinear to initial-state partons, 



rl *N-\ _ i T r (l-z)Q 2 Arp. 

\nJ q (Q 2 ) = \ dz— / J^ A (a s (q 2 )) + B(a s ([l-z]Q 2 )) 

Jo 1 — 2 J(i-z) 2 Q 2 q z 



(4) 



<±-z)Q 2 dq 
l(i-z) 2 Q 2 q 

taking into account collinear final-state radiation, and the process-dependent piece 

r i z n-i _ i 

lnA^(Q 2 ) = / dz— D P (a s ([l-z] 2 Q 2 )) (5) 

Jo 1 — z 

attributed to large-angle soft-gluon emissions. The integrands in eqs. ©-(0) are given by 

oo 

F(a s ) = F = A,B,D P , (6) 



i=i 



1 



where we normalize the expansion parameter as a s = a s j (Air) . 

The constants A\ in eq. @ are the coefficients of the 1/[1 — x] + terms of /-loop quark- 
quark splitting functions P^~ x \x) — this actually completes the MS definition of the 
coefficients Bi and Dpi (see also ref. ||). Thus A% and A 2 are well-known, reading 

'67 u \ „ 5 
9 



A l = AC F , A 2 = 8C f 



C 9 )C A --N f 



(7) 



Here Nf denotes the number of effectively massless quark flavours, and the colour factors 
are Cf = 4/3 and Ca = 3 in QCD. The exact expression for P@'(x) has not been 
completed yet. However, recently rather accurate approximations have been derived |J 
from the available partial results, most notably the lowest integer- iV moments [[|. |8|]. 
Together with the exact N 2 term determined in ref. @ the results of ref. || yielcf] 

16 

A 3 = (1178.8 ± 11.5) - (183.95 ± 0.85) N f -— C f N 2 f . (8) 

27 3 

For Nf = 3 ... 5 the (independent) errors in eq. (|§|) can be combined to an overall 
uncertainty of ±12. For a s < 0.3 this uncertainty amounts to less than 0.1% of the total 
three-loop value of A(a s ). Finally the constants E>i and Dpi in eq. @ are given by 

Bi = -ifi = -3C F , Dpa = , (9) 

where P j is the coefficient of 8(1 — x) in the one- loop quark-quark splitting function. 
The second-order terms B 2 and D P2 will be discussed below. 

After the integrations in eqs. (ffH© are performed, the functions Gp{Q 2 ) in eq. (fj) 
take the form 

G N (Q 2 ) = Lg 1 (X)+g 2 (X)+a s g 3 (X) + ... (10) 
with L = liaN, X = /3oa s L and 

00 

9iW = HdAasLf . (11) 

k=l 

The first term in eq. (PD, which depend on A\ only (see eq. (|TB|) below), collects the leading 
logarithmic (LL) large-iV contributions L(a s L) n . The coefficients A 2 and B\ determine 
the functions g 2 resumming the next-to-leading logarithmic (NLL) terms (a s L) n . These 
functions have been determined in refs. [|1], The next-to-next-to- leading logarithmic 
(NNLL) approximation includes g 3 which is correspondingly fixed by A 3 , B 2 and Dp 2 . 

In order to calculate g 3 it is convenient to introduce 

a 2 

X = l + a a (i4)p hi±= (12) 

Mr 



■'^The error estimate m) has been derived in ref. under an assumption (no ln 3 (l — x) terms) on the 
form of Pqq (x). Abandoning this constraint does not lead to a significant modification of eq. (H). 
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and to write the next-to-next-to-leading order coupling constant as 



X 



2 2 fa lnX 



3 2 fPfln 2 X -\nX -1+X fa 1 - X 



+ o>:(j4) ik 



+ ~77~ 



(13) 



M X3 ' 0o X* J 

+ 0(a 4 s (fil)[a s (fi 2 r )\n(q 2 /^ r )] n ) . 

Here /3q, Pi and Pi are the coefficients of the /3-function of QCD up to three loops jlOf . 
After inserting eq. (|13"D into eqs. @ and (f|), the inner integrations over A(a s ) can be 
readily carried out. A straightforward method to perform the ^-integral is to write the 
integrand as an infinite series in ln(l — x), to use 



JV-l 



dz 



1-z 



ln fc (l - x) 



[ -^- {^ +1 (iV) + \ k(k + 1) Sf\N) S 2 (N) 



+ o^- 2 ) 

(cf. ref. 0|) with S a (N) = EjLi l/j a and 

^(iV) = lniV + 7e + 0(l/iV), S 2 (N) = ( 2 + (D(l/N) ., 
and to re-assemble the expansions in the end. In this way we arrive at 
nis/x\ 



9 rw 

# S (A) 



p x 



A+(l- A)ln(l-A) 



(14) 



(15) 



(16) 



^e-Bi_ Hl _ x) + AiPi 



01 



Po 

A + ln(l - A) 



A + ln(l-A) + ^ln 2 (l-A) 



+ In 



Mr 



and our new result 
^ 3 DIS (A) 



, I 1 f 2 ^ X 01 1 



1 



lln 2 (l-A) + Aln(l-A) + -A 2 



01 1-A 
A 



A + ln(l - A) 



+ 



/?o 3 



1 A 2 



A 1 
Po 1 - A ~ /3q 1 - A 



2 1 - A 
1 



ln(l-A) + A + ^A 2 



+ ln(l-A) + A 
> + * A2 



- Bi 7e 



A 



ft 



A p 2 1 - A 



A + ln(l - A) 



2/3 2 1 - A 
-B2 A D12 A 



A) 1 - A fa 1 - 2A 



/3 2 1 -A 



Q 2 \ A 



+ W -i 



A*? 



2 1-A 



A + ma-A^^-^-f)^ 



H 2 - J Po 



Mr 



(18) 



The functions g 2 J are obtained from eqs. ( |T7D and (|TJ) by substituting A — > 2A, 7 e — > 27 e 
and — » 4(^2 in the terms with Aj, removing the terms with (recall eq. fl2|)), and 
replacing D11 by D 2) i. The result corresponding to eq. (P3|) reads gf Y (\) = 2gf IS (2A). 
The generalization of eqs. (|16D — (|T~8D to other processes involving eqs. ©-(§D is obvious. 

Now we are ready to address the second-order coefficients B 2 and Dp 2 . The first-order 
expansion coefficients gfP and gf Y , as defined in eq. (|TT|), are given by 



(P) 

031 



1/2 P 3 A 1 ( 7e 2 + C 2 )/?o + P 2 A 2le - 6 P1 (B 2 - B lle (3 ) - D P;2 



(19) 



(p) 

with 5kj = 1 for k = j and 8kj = else. On the other hand can be determined 
by expanding eqs. (P and (jlOf) to order a 2 and comparing to the two-loop results of 
refs. |3|, |T2|] , as done for the DIS case in ref. ||1_3[Q The result for the Drell-Yan coefficient 
function reads 



9^ 



-56C 3 
32 2 

T % 2 + 



176 „ 1072 

32C 2 7e + — 7e 2 + -^- 



7t 



160 



yielding 



DY 



CfCa 



1616 r „. 176 . 



7< 



/224 32 
+ C F N f (-- Y Cz 



(20) 



(21) 



This term has already been derived in ref. [14], albeit without explicitly attributing it to 
the a s ((l — z) 2 Q 2 ) contribution. Note that, as it has to be for a A int contribution M, 
eq. ( pi] ) does not contain an Abelian (Cp) term. We consider this as a first concrete check 
of the correctness of the setup at the NNLL level. Using eq. (12) of ref. |L3|, the 

corresponding constraint for the DIS case reads 



B 2 + D 2 DIS 



where P^j is the coefficient of 8(1 — x) in the two- loop quark-quark splitting function. 
Unlike the DY case fl21~|) two new constants can occur at the NNLL level in DIS, hence the 
consistency of the resummed and the two-loop coefficient functions does not completely 
specify g^ ls - For that either an extension of the calculations of refs. |1], |2| to the next 
order, or the In 2 iV term of the three-loop coefficient function (fixing g^ 2 which involves 
the combination B 2 + 2_D lj2 ) is required. An approximate result has been derived for the 
latter [15 from the constraints of refs. |7], || 13fl . Within errors that result is consistent 



--24C3 + 12C2 



/247 
+ C F N f (- 



+ c F c A - 



3155 
~54~ 



1 
2 



44 \ 
+ 40C 3 + yC 2 J 

7 (3 Cf , 



'.22) 



with D® 1S 



0, but, for instance, also with Bo 



(i) 



q,S 



£p C F for £ ~ 8 ... 13. 

2 Note that for 173 the convention differs from that in ref. [|l3|: there the second index i in g$i 
refers to the overall power of a s in the expansion (jl0|) of G N . Hence 1731 is called 332 in ref. 1 1 3 etc. 
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The integrals in eqs. (|3|)-(|5|) are not well-defined, as they involve the running coupling 
at arbitrarily low scales. This feature leads to the poles at (2) a s /5 ln N = 1 in eqs. (fL6f)- 
and their DY counterparts, a problem which is usually dealt with by the 'minimal- 
prescription' contour for the Mellin inversion [16]. Another option, followed in ref. jDJ 
for the DIS case, is to re-expand eq. ([!]) in powers of a s and to keep only those terms 
a™ \n 2n+1 ~ k jV, k = I,... I (i.e., the first I 'towers' of logarithms) which are completely 
fixed by the known terms in eq. (10) and in go of eq. (|l|). In connection with a two- loop 
result for go the NLL and NNLL resummations lead to / = 4 and / = 5, respectively. 
General expressions for the first four towers can be found in eq. (14) of ref. |13| 2 , the 
extension to higher terms is straightforward if lengthy. After Mellin inversion the NNLL 
resummation thus leads to the following prediction for the five leading x — > 1 terms 
of the three-loop (quark-antiquark annihilation) DY coefficient functions for dojdQ 2 at 
fij = fif, = Q 2 (the other terms are fixed by renormalizat ion-group constraints): 



J3) 



512 C% 



ln 5 n -x) 



(7 _m _ mo clN 



V 9 9 
/17152 



p f 



ln 4 (l - x) 



+ {(- 2048 - 3072 C 2 ) C% + (—^- - 512 C 2 ) C 2 C A + 



1 — x 
7744 



2560 ^ 2 Ar 2816 _ _ Ar 256 „ Ar2 

ClNr C F C A N f + C F N 2 f 

9 F f 27 F A f 27 F f 



27 

ln 3 (l -x) 



C F C\ 



1 — x 



+ { 10240 C 3 G% + ( - ^ + 1344 C 3 + ^C 2 ) C\C A 



28480 . 704 
3 



+ 



27 

/9248 128 

\in 3~ 



' )c ^ ^644 2048 , lrf 



V 9 

C 2 )C F C A N f -^C F N 2 f 



ln 2 (l -x) 
1 — x 



(23) 



-( 



77949.50 - 5886.63 N 



/ 



32.888 Nf + AA 3 



ln(l — x) 
1 — x 



+ 



Here c^ Y (x) refers to a s = a s /(47r), as all other expansion parameters in this paper. 
Since A 3 is presently known just in the form (§), only the numerical result has been 
written down for the [(1 — a;) _1 ln(l — x)] + term in eq. (|2"3"D. In practice, however, the 
uncertainty of ±12 for A% is irrelevant here and in the corresponding higher-order terms. 



B 2 1S in eq 



An improvement on the tower expansion for DIS in ref. will only be possible once 
^2|) is exactly determined. However, the present results for g^ ls are directly 
relevant to the 'physical' evolution kernels K N for the scaling violations of the non-singlet 
structure functions investigated beyond order a 2 in refs. |T5|, [H]]. Identifying all scales for 



3 Note that the third term of eq. (15) has been misprinted in the original preprint as well as the journal 
version of ref . @: 16/27 C F N f has to be replaced by W/27C F Nj, and 70/9C F Nj by 280/9 CjJV/. 
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brevity these kernels can be written as 

= (P N (a s ) + dlnC d : (as) P(*s)) F a %(Q 2 ) = K*(a s ) F a %(Q 2 ) , (24) 



where P N (a s ) and C^(a s ) (a = 1, 2, 3) are the moments of the non-singlet splitting 
functions and coefficient functions, respectively. At large-x / large- N eq. fl2"4|) holds for 
the full structure functions up to corrections of order 1/N. Inserting eqs. (|1|) and ( [T0| ) 
into this equation one arrives at the NNLL resummation of this kernel, 

KM = -\nN(A ias + A 2 al + A 3 al) - (l + ^a s + ^\ 2 ^ (25) 



A + a%)\^ - a 2 s (3 -^-(Xg 3 (X)) + 0(a 3 s (a s \nNy 



dX dX 

with A = a s P \nN, A4 of eqs. (|7|) and (||), and g i from eqs. (|Ilf)-(|i"8|) with /// = \i T = Q. 
Note that the leading large- N term of the quark-quark splitting function is given by In N 
at all orders in a s [|T8]]. 

Finally we briefly illustrate the numerical impact of the NNLL corrections. For these 
illustrations we choose fi 2 = fi 2 = Q 2 , Nf = 4 and a s = 0.2. Depending on the precise 
value of a s (M|), the latter number corresponds to scales between about 25 and 50 GeV 2 , 
a range typical for fixed-target experiments both on DIS and the Drell-Yan process. The 
corresponding results for the functions G^ 1S and G^ Y are presented in fig. 1. Their 
convolutions with a schematic, but typical input are shown in fig. 2, where the minimal- 
prescription Mellin inversion |Tj| has been employed. The remaining ambiguity in the 
DIS case (discussed below eq. (|2|)) is illustrated by the results for -DJ? IS = and for 
B 2 = -P$ + 13 (3 C F (denoted by f = 13 in the fi gures), the latter curves indicating the 
maximal NNLL effects consistent with the available three-loop information [[F], |], |H| [15 



It is obvious from both figures that knowledge of the In N(a s In N) n and (a s In iV) n terms 
[l], alone is not sufficient for reliably determining the functions Gp and their impact 
after convolution even for rather moderate values of iV and x (often denoted r in the 
DY case). On the other hand the NNLL corrections presented in this paper are rather 
small over a wide range, e.g., less than 10% and 20% at x < 0.85 and y/x < 0.75 for the 
DIS and DY results of fig. 2, respectively. This stabilization indicates that the soft-gluon 
exponents Gp and their effects can now be reliably estimated for these processes. 
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Fig. 1: The LL, NLL and NNLL approximations for the resummation exponents G N (Q 2 ) 
in eq. (|H]) at p? r = \A = Q 2 for a s {Q 2 ) = 0.2 and four flavours. The two NNLL curves in 
the DIS case indicate the present uncertainty as discussed below eq. fl22|). 
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Fig. 2: The convolutions of the results shown in Fig. 1 with a typical input shape. 
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